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Abstract 

In this paper we study fj, — e conversion in nuclei within the context of the Con- 
strained Minimal Supersymmetric Standard Model, enlarged by three right handed 
neutrinos and their supersymmetric partners, and where the neutrino masses are gen- 
erated via a seesaw mechanism. Two different scenarios with either universal or non- 
universal soft supersymmetry breaking Higgs masses at the gauge coupling unification 
scale are considered. In the first part we present a complete one-loop computation of 
the conversion rate for this process that includes the photon-, Z-boson, and Higgs- 
boson penguins, as well as box diagrams, and compare their size in the two considered 
scenarios. Then, in these two scenarios we analyse the relevance of the various param- 
eters on the conversion rates, particularly emphasising the role played by the heavy 
neutrino masses, tan/3, and especially ^13. In the case of hierachical heavy neutrinos, 
an extremely high sensitivity of the rates to ^13 is indeed found. The last part of this 
work is devoted to the study of the interesting loss of correlation between the /i — e 
conversion and ji ^ rates that occurs in the non-universal scenario. In the case of 
large tan (5 and light Higgs boson, an enhanced ratio of the ;U — e to — > 67 rates, 
with respect to the universal case is found, and this could be tested with the future 
experimental sensitivities. 



1 Introduction 



Neutrino physics, in particular neutrino oscillations and the measured neutrino mass differ- 
ences, strongly manifest that Nature does not conserve the lepton flavour quantum number 
in the neutrino sector [1, 2]. However, it is not known yet if lepton flavour violation (LFV) 
also occurs in the charged lepton sector. If such is the case, one still has to address if LFV 
in the neutral and charged lepton sectors arises from a common or different origin. It is well 
known that if the Standard Model of Particle Physics (SM) is minimally extended in order 
to accommodate the present data on neutrino masses and mixings, the corresponding loop 
induced LFV in the charged lepton sector (exclusively induced from neutrino oscillations) 
is extremely tiny and hopeless to be experimentally observed. Therefore, a potential future 
measurement of LFV in the charged lepton sector will provide a unique insight into the 
nature of new physics beyond the SM (for a review, see [3]). 

Among the various candidates for physics beyond the SM that produce potentially ob- 
servable effects in LFV processes, one of the most appealing are Supersymmetric (SUSY) 
extensions of the SM, where a seesaw mechanism [4, 5] is implemented to generate neutrino 
masses. In these SUSY-seesaw models a new source of LFV appears in the off-diagonal 
elements of the slepton mass matrices, which can be radiatively generated. The size of these 
elements is governed by the strength of the neutrino Yukawa couplings and, in the case of 
Majorana neutrinos, the latter can be large, of the order of one. The LFV effects in the 
charged lepton processes arc then induced by flavour violating slepton-lepton interactions, 
appearing in SUSY-loop diagrams mediated by sleptons [6]. 

Concerning the LFV processes, in our work we will focus on those which involve flavour 
transitions between the first and second generation of charged leptons. At present, the most 
relevant — e flavour violating processes are /j, 67, /j, ^ 3e and ^ — e conversion in nuclei. 
The current experimental bounds on the muon decays are BR(// 67) < 1.2 x 10^^^ [7] and 
BR(yU 3e) < 1.0 X 10^^^ [8]. Regarding fi — e conversion in heavy nuclei, the most stringent 
constraints arise for Titanium and Gold, respectively with CR(/i — e, Ti)< 4.3 x 10^^^ [9] and 
CR(/i — e, Au)< 7 X 10^^'^ [10]. In the future, one expects significant improvements in the 
sentitivies to these LFV rates. For instance, MEG aims at reaching a sensitivity for fi ^ e'-f 
of 10~^^ [11] in the very near future, which could further be improved to 10~^^ in the next 
4-5 years [12]. Although the situation for BR(/x 3e) is less certain, one does not expect 
the sensitivities to better 10^"*^'^ — 10^"*^^ [12]. Undoubtedly, the most challenging prospects 
concern the experimental sensitivities to fi — e conversion in Titanium nuclei. The dedicated 
J-PARC experiment PRISM/PRIME has anounced a remarkable improvement, albeit in a 
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farer future, of 10"^^ [13]. 

In this paper we will focus on // — e conversion in nuclei, working in the context of the 
Minimal Super symmetric Standard Model (MSSM) enlarged by three right handed neutrinos 
and their corresponding SUSY partners, where a type-I seesaw mechanism [4] is implemented. 
To reduce the number of unknown parameters in the SUSY sector, we choose to work in 
the so-called constrained MSSM (CMSSM)(for a review see for instance [14]), assuming 
universahty of the soft-SUSY breaking parameters at the scale of gauge coupling unification, 
Mx ~ 2 X 10^^ GeV. An interesting departure from the CMSSM-seesaw can be obtained 
by relaxing the universality hypothesis for the soft SUSY breaking masses of the Higgs 
sector. This partially constrained MSSM is commonly referred to as the Non Universal Higgs 
Mass (NUHM) scenario [15], and its enlarged version (including right handed neutrinos and 
sneutrinos) will be here designated NUHM-seesaw. 

Within the context of the CMSSM- and NUHM-seesaw, we conduct here a thourough 
analysis of the predictions for the /i — e conversion rates in nuclei. The present computation 
is the first to include the full set of SUSY one-loop diagrams (photon, Z- and Higgs-boson 
mediated, as well as box diagrams), and to be strictly done in terms of physical eigenstates 
for the exchanged SUSY particles. For all scenarios here addressed, we obtain the low-energy 
parameters by numerically solving the full renormalisation group equations (RGEs), includ- 
ing the neutrino and sneutrino sectors. The photon and Z-boson mediated penguins, as well 
as the vector contributions from box diagrams were first computed in the CMSSM-seesaw 
in [16]. Here we have confirmed their analytical results for the photon- mediated and box 
contributions, correcting the analytical expressions for the Z-boson mediated processes. We 
further added the scalar contributions from box diagrams and the Higgs-mediated contribu- 
tions, and improved the computation, by considering in the numerical analysis the possibility 
of either degenerate or hierarchical heavy neutrino spectrum, and by fitting the light neutrino 
parameters to the present data. 

The effects of the Higgs-mediated contribution on // — e conversion rates were firstly 
investigated in [17], in the context of a SUSY-seesaw with degenerate heavy neutrinos, 
working in the effective Lagrangian approximation and in the SU{2)l x U{1)y limit. It was 
observed that the Higgs-mediated LFV diagrams could provide the dominant contribution 
for the large tan /3 regime and for small masses of the heavy Higgs scalar, owing to a tan^ (3 
enhancement and (m/fo)"^ dependence of the conversion rates. By comparing the latter with 
the corresponding — > 67 rates for large universal SUSY-breaking mass values, Mq, M1/2 ~ 
0{1 TeV), they further showed that the ratio of observables CR(/i. — e, Al)/BR(yU —>■ ej) 
could be enhanced from a value of 0{a) (within the usual dominant photon- mediated case) 
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to 0(10-^) for extreme values of tan/3 = 60, Mr = 10^^ GeV and m^jo ~ 100 GeV. 

In the present work we will explore in full detail the various contributions to the — e 
conversion rates and study the dependence on all parameters entering in the considered 
MSSM-seesaw framework. In addition to the relevant role played by the mass of the right 
handed neutrinos, mAr^, the soft masses Mq, M1/2 (and M^^ j for the NUHM-seesaw) , and 
tan/3, we will show that the light neutrino mixing angle ^13 has an important impact on 
CR(/x — e. Nucleus). The conversion rates turn out to be very sensitive to this angle, varying 
in many orders of magnitude (up to five in the case of hierarchical heavy neutrinos) for ^13 
values within the present experimentally allowed region 0° < ^13 < 10° [18]. We will further 
verify that with the future sensitivity of JPARC (O(10~^^)) [13] most of the parameter space 
could be covered. 

On the other hand, the comparison between the predictions obtained for the CMSSM- 
seesaw and the NUHM-seesaw cases will allow us to draw interesting conclusions about the 
departure from the strongly correlated behaviour of CR(// — e. Nucleus) and BR(// — > 67), 
as predicted in photon-dominated scenarios (as is the case of the CMSSM-seesaw) . In the 
latter scenario, the ratio of the two rates was found to be at most 1/160 for tan /3 = 50 [19]. 
In contrast, we will discuss here particular scenarios in the NUHM-seesaw, where the ratio 
CR(/x — e, Ti) /BR(// — > 67) is indeed enhanced with respect to the universal case, by as 
much as one order of magnitude, in agreement with the approximate results of [17]. 

One of the most challenging tasks in this ^ — e conversion process will be to disantangle 
between the different scenarios for new physics if a measurement is finally obtained. Indeed, 
it has been already noticed in early works [20] that ^ — e conversion could constrain new 
physics more stringently than ji — > 67. We will see here that this is the case in the NUHM 
scenario. Furthermore, we will also show that, with the expected sensitivities for Titanium 
of 0(10"-*^^), one could distinguish CMSSM- from NUHM- seesaw scenarios by extracting 
the scalar contribution to the CR(// — e, Ti) rates. 

The paper is organised as follows. In Section 2 we review the most relevant features of the 
SUSY-seesaw scenario. The analytical results of the ji — e conversion rates are presented in 
Section 3. Section 4 is devoted to the numerical results for both CMSSM-seesaw and NUHM- 
seesaw scenarios. An extensive discussion about the sensitivities to the various parameters 
in these two scenarios is also included. Finally, the conclussions are summarised in Section 5. 



4 



2 The SUSY-seesaw scenario 



The leptonic superpotential containing the relevant terms to describe a type-I SUSY seesaw- 
is given by 

W = N^Y^LH2 + E^YiLHi + ^N^niMN^ (1) 

where arc the additional supcrficlds that contain the three right-handed neutrinos v^^ 
and their scalar partners v^.. The lepton Yukawa couplings Yi^y and the Majorana mass mu 
are 3x3 matrices in lepton flavour space. From now on, we will assume that we are in a 
basis where Yj and are diagonal. 

After clcctroweak (EW) symmetry breaking, the charged lepton and Dirac neutrino mass 
matrices can be written as 

= f 1 , mo ^Y^V2, (2) 

where Vi are the vacuum expectation values (VEVs) of the neutral Higgs scalars, with vi(2) = 
V cos(sin)/? and v — 174 GeV. 

The 6x6 neutrino mass matrix is given by 



m 



M'^ = ^ . (3) 

rriB vhm 

The eigenvalues of are the masses of the six physical Majorana neutrinos. In the seesaw 
limit, the three right-handed masses are much heavier than the EW scale, ^ v, and 
one obtains three light and three heavy states, and N^, respectively. 

Block-diagonalisation of the neutrino mass matrix of Eq. (3), leads (at lowest order in the 
{mB/mMY expansion) to the standard seesaw equation for the light neutrino mass matrix, 

= -m^rn^mD ■ (4) 

Since we are working in a basis where tum is diagonal, the heavy eigenstates are then given 
by 

m%^^ = tum = diag (mj^^, mN2, ruN-i) ■ (5) 

The matrix rriiy can be diagonalised by the Maki-Nakagawa-Sakata unitary matrix C/mns [21, 
22], leading to the following masses for the light physical states 

mt'^^ = f^MNS "^i- ^MNS = diag{m^^,m^^,m^^) . (6) 
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Here we use the standard parameterisation for C/mns given by 



Umns 



Cl2 Ci3 
-Sl2 C23 - C12 S23 -§13 e*' 



S12 Ci3 
C12C23 - Sl2S23Sl3e* 



V S12S23 - C12 C23 si3 e*'' -C12S23 - S12 C23 si3 e*'' 



si3 e-'^ \ 

S23 Ci3 
C23 Cl3 / 



with 



V 



V. (7) 



(8) 



diag (e ' 2 , e * ^ , 1) , 

where = cos 9ij, Sij = sin6'ij. % are the neutrino flavour mixing angles, 5 is the Dirac 
phase and ^1^2 are the Majorana phases. 

In view of the above, the seesaw equation (4) can be solved for mD as [23] 



mo 



i\lmf^R\lmf''^U} 



MNS ' 



(9) 



where i? is a generic complex orthogonal 3x3 matrix that encodes the possible extra neutrino 
mixings (associated with the right-handed sector) in addition to the ones in C/mns- can 
be parameterised in terms of three complex angles, (i = 1, 2, 3) as [23] 



R 



^ C2 C3 -Ci S3 - Si S2 C3 Si S3 - Ci S2 C3 

C2 S3 Ci C3 - Si S2 S3 -Si C3 - Ci S2 S3 
\ "52 Si C2 Ci C2 



(10) 



with Cj = cos^j. Si = sin^j. Eq. (9) is a convenient means of parameterising our ignorance 
of the full neutrino Yukawa couphngs, while at the same time allowing to accommodate 
the experimental data. Notice that it is only valid at the right-handed neutrino scales 
niM, so that the quantities appearing in Eq. (9) are the renormalised ones, m^'^^ {ruM) and 
C^MNS (itim)- 

We shall focus on the scenario where the light neutrinos are hierarchical, and we will 
assume a normal hierarchy. 



(11) 



The masses m,/2,3 can be written in terms of the lightest mass m.i,^, and of the solar and 
atmospheric mass-squared differences as 



m 



1/2 

m^3 = Am^tm + 



(12) 
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Regarding the heavy neutrinos, we will consider the two following cases, 

Degenerate: m^r^ = 171^^ = rrij^^ =171^, 
Hierarchical: ttini ^ '>tin2 ^ miVg ■ 

Concerning the SUSY parameters, and since we are working within an extended MSSM, 
with enlarged neutrino and sneutrino sectors, there will be new soft SUSY breaking parame- 
ters associated to the latter sectors. Thus, in addition to the usual soft breaking parameters 
for the gauginos (Mi^2,3), Higgs bosons (M//^ ^); squarks {ttiq, rrifj, mi), Ag) and sleptons {mi, 
mi;, Ai), there will also be the sneutrino soft breaking masses m^g, the sneutrino trilinear 
couplings A^, and the new bilinear parameter Bm- As already mentioned in the introduction, 
we will work in a constrained MSSM, where the number of input parameters is reduced by 
assuming partial universality of the soft parameters at the gauge coupling unification scale, 
Mx — 2 X 10^^ GeV. Specifically, we will work in two scenarios, the CMSSM-seesaw with 
universal scalar masses, trilinear couplings and gaugino masses, and the NUHM-seesaw with 
non-universal soft masses for the Higgs bosons. Therefore, when specifying the parameters 
of these two constrained MSSM scenarios we will fix, in addition to the seesaw parameters, 
the following soft SUSY breaking parameters at the scale Mx- 

CMSSM-seesaw: Mq , M1/2 , , tan /5 , sign(/i) , 
NUHM-seesaw: Mq , M1/2 , ^0 , tan /? , sign(//) , Mh^ , Mh^ . (13) 

The departure from universality in the NUHM-seesaw will be parameterised in terms of the 
non- vanishing parameters 5i and 82, 

Non-universality: = (1 + 81) , Mjj^ = Ml (1 + ^2) • (14) 

For simplicity, and to further reduce the number of input parameters, in this case we will 
also impose Mq = M1/2 = -^susy- 

Once the above set of parameters is fixed at Mx, the predictions for the low-energy 
parameters are obtained by solving the full one-loop RGEs, including the extended neutrino 
and sneutrino sectors. Due to the existence of intermediate scales mu introduced by the 
seesaw mechanism, the running must be carried out in two steps. The full set of equations 
is first run down from Mx to mM- At the seesaw scales, the right-handed neutrinos as well 
as their SUSY partners decouple, and the new RGEs (without the equations and terms for 
vr and vr) are then run down from mu to the EW scale, where the couplings and mass 
matrices are finally computed. 
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Working in constrained MSSM scenarios, all flavour mixing originates solely from the 
neutrino Yukawa couplings, which induce flavour violation in the slepton sector by the RGE 
running from Mx down to the EW scale niz- Flavour mixing is then manifest in the values 
of the off-diagonal elements of the charged slepton squared mass matrix. The LL, RR, LR 
and RL elements of the latter matrix can be summarised as follows: 

M^L = rnl .. + v\ + m| cos 2/3 (^-^ + sin^^i^^ 5,,-, 

M'il = + v\ - ml cos 2/3 sin^^^y^,,. 

Mil = (Mf^^)* . (15) 

In the above, mz denotes the Z-boson mass, 6w the weak mixing angle, and = 1, 2, 3 are 
flavour indices. Given that below itlm the right-handed sncutrinos decouple, the low-energy 
sneutrino mass eigenstates are dominated by the vl components. Thus, sneutrino flavour 
mixing is confined to the left-handed sector, and described by the following 3x3 matrix: 

Mr = m|_^^. + ^m|cos2/35,,. (16) 

The physical masses and states are obtained by diagonalising the previous mass matrices, 
leading to 

^2diag ^ ^(/)^2^(0t = diag(m?,..,m?), 

^2diag ^ ^(.) ^2 ^(.) t ^ diag {ml , ml , , (17) 

where R^^''^^ are unitary rotation matrices. 

After having introduced our scenario, in the next section we will summarise some of the 
more relevant details leading to the computation oi fi — e conversion rates in nuclei. 



3 Analytical results of the ji — e conversion rates 

In this section we report the analytical results for the — e conversion rates in terms of the 
parameters introduced in Section 2. We emphasise that all the results are obtained in terms 
of physical mass eigenstates (with full propagators) for all MSSM particles entering in the 
computation, namely, charginos xli-^ — 1; 2), neutralinos x^(A = 1, 4), charged sleptons 
lx{X — 1, 6), sneutrinos Vxi^ — 1> 2, 3) and the neutral Higgs bosons, hP and 
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For the presentation of the results we closely follow the general parameterisation (and 
approximations) of [3]. One starts with the most general effective Lagrangian for four- 
fermion interactions which describes coherent — e conversion. At the quark level, this is 
given by 

Gp . 

-^eff = -J= ^ { [gLS{q)eLlJ^R + 9RS{,q)eRHL] QQ + [q LV {q)^ L + 9 RV {q)^ R-^^ ^ r\ QlfiQ} , 

q 

(18) 

where is the Fermi coupling. Notice that only scalar (S) and vector (V) effective operators 
do contribute, with couphngs given by gLS{q), gRS{q) and gLV{q), QRViq) (respectively left and 
right, in both cases). This effective Lagrangian at the quark level is then converted into 
an effective Lagrangian at the nuclcon level, by means of the appropriate nucleon form 
factors [24]. In the limit of negligible momentum dependence of the nuclcon form factors, 
(a reasonable approximation given the small momentum transfer in the fi — e process), the 
quark matrix elements can be simply replaced by the nucleon matrix elements as follows: 

{n\qrKq\n)^G^r^nrKn, (19) 

where Tk — (1,7/*) respectively for K — {S, V). The numerical values of the relevant Gk^s 
are [3,25]: 

q{u,p) _ Q{d,n) _ 5 . Q{d,p) _ Q{u,n) _ ^ g . 

S S ' ' iS S ' ' 

q{s,p) ^ ^is,n) ^ . (20) 
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The conversion rates are then predicted in terms of the relevant isoscalar, g^xK^ isovector 
couplings, g^xK (with X — L,R and K — S,V), which are given by: 

q=u,d,s 

& - \ E {9xK^.)G^r^ - gxKiq)G^r^) . (21) 

q=u,d,s 

Further working under the approximation of equal proton and neutron densities in the 
nucleus, and of a non-relativistic muon wave function for the 1 s state, the final formula for 
the /X — e conversion rate, relative to the the muon capture rate, can be finally written as 

CR /x - e, Nucleus = " 

x<(Z + N) (ji»> + +(Z-N) + «S) f + 

(Z + N) (sS + sis) +(Z-N) (sS + siS) f } . (22) 

where Z and N are the number of protons and neutrons in the nucleus, while Z^s is an 
effective atomic charge, obtained by averaging the muon wave function over the nuclear 
density [26]. Fp is the nuclear matrix element and Fcapt denotes the total muon capture rate. 
The other quantities in the above formula correspond to the muon mass, m^, the momentum 
and energy of the electron, pg and (which are set to in the numerical evaluation), and 
the electromagnetic coupling constant, a. 

We have computed the full set of one-loop diagrams leading into the quantity CR(/i — 
e. Nucleus): 7-penguins, Z- and Higgs-boson penguins and box diagrams. These are schemat- 
ically drawn at the quark level in Fig. 1, and receive contributions from several diagrams, 
mediated by SUSY particles, which are collected in Appendix A. The analytical results 
of the computation are summarised in terms of the contributions of these diagrams to the 
vector and scalar couplings, 

gLV{q) = gLV{q) + 9LV{q) + gLV{q) > 

gLS{q) = gLS{q) + 9lv{q) " (23) 

In the above, the photon couplings g1x{q)^ Z-boson couplings gfx(q)' i^-boson cou- 
plings g^s{q)i couplings arising from the boxes g^x{q) (with X = V,S) are respectively 
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given by 



9LV{q) ~ 'G^^^^ ^2^^ ' 

9LVi,)- 2 

_ fr>{n)LV , r>ic)LV\ 

9LV{q) - -JT +^q' ) 

V2I 



M _ V2 1 ^-^ 1 (p) / (p) 

^ p Hp 



= + (24) 

Likewise, for the right-handed couphngs we find 

9RV{q) = 9LV{g)\L^ji , 

gnsiq) = 9LSiq)\L^ji ■ (25) 

The exphcit formulae for the form factors of the photon {A^^'^^ ), of the Z-boson (-F(l,_r)), 
of the Higgs-boson {Hj^^j^y where p = 1,2,3 corresponds to Hp = h^^H^^A^)^ and of the 
box diagrams (^b'^'^^'^^'^'^^^'^'') are hsted in Appendix A. In each case, the relevant couplings 
^{LR)^ ^'^R)q ^® found in Appendix B. 

It is important to stress that 5"^^^ + 5"^^^ vanishes and therefore there are no contributions 
from the CP-odd Higgs boson A^. This is a consequence of working in the approximation of 
coherent // — e conversion, in which case the initial and final nucleus state is the same, thus 
leading to vanishing matrix elements for pseudoscalar currents like (Nucleus | g 75 g | Nucleus) . 
Also notice that from the values of the S'^j^-^^ Higgs couphngs, one can anticipate that in 
the large tan /3 and small Higgs mass regime, the dominant Higgs contribution will be that 
of H°. 

When compared to the results obtained in [16], our expressions coincide in the formulae 
for the photon-penguins. Up to a global sign, the vector contributions from boxes also 
agree. Divergences occur regarding the Z-penguins, and the differences can be read by 
comparing our expressions in Eqs. (44-48) of Appendix A, with those of Eqs. (22-29) in [16]. 
As previously mentioned, we have included in addition scalar contributions from boxes and 
Higgs- mediated diagrams not considered in [16]. 

A connection between our results for the Higgs contributions and those reported in [17] 
can be established in the large tan /3 limit, writing the output in the mass-insertion approx- 
imation format. Under these conditions, and considering the limit of a common mass for all 
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SUSY particles involved, which is much larger than the SM particle masses, Mgoft ~ A* 
-^susY >> 'iTT'W, one arrives at the following simple expression for the dominant form 
factor [27] 



^ (47r)2 12mvK 



1 + -(1 -Stan^ ^w^) 
2 



(26) 



where the first term arises from chargino mediated loops, while the second stems from neu- 
tralino mediated contributions. In the mass insertion approximation, the dominant slepton 
mixing effects are associated with ^21) which can be written as: 

(Am?) 21 

-'^^SUSY 

From the above, one can finally obtain a simple expression for the contribution to the 
conversion rate, which is clearly dominated by the strange quark coupling, due to the en- 
hancement in the coupling by rus- This arises via g^^l ~ 9ls{s)^s'^^ with 

HO _ V2 1 1 (2) gms 
^^^(^) -G^2^,^^ ^^^^ 

Plugging this simplified result for the g^^l coupling into the approximate conversion rate for 
the Higgs-dominated case. 



CR(/. - e. Nucleus) "g^^^ ' iiZ + N) g^' \ ^ (29) 



Tcapt 



we obtain the expected tan^ /3 enhancement of the contribution. Moreover, the depen- 
dence on the Higgs mass i:^^)-, as well as the typical pref actor l^ail^ accounting for the 
lepton flavour changing effect are equally recovered. Within this approximation, and tak- 
ing a specific value of 5\-^ — 10~^, allows to obtain an order-of-magnitude estimate for the 
conversion rate in the case of Titanium nuclei, 

115GcV\^ /tan/? 



in agreement with the estimate of [17]. 

Finally, it is worth mentioning that the heavy SUSY particles do not decouple in the Higgs 
contributions to the — e conversion rates. This can be understood from the previous result 
of in Eq. (26), which is constant in the large Msusy limit. This SUSY non-decoupling 
effect has also been noticed in association to other Higgs- mediated LFV processes [27-30]. 
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4 Numerical results and discussion 



In this section we present the numerical results of the — e conversion rates in nuclei within 
the SUSY-seesaw context described in Section 2. We begin by addressing the CMSSM- 
seesaw, and then proceed to the NUHM-seesaw. In both scenarios, we consider the depen- 
dence of the theoretical predictions for the conversion rates on the most relevant SUSY- 
seesaw parameters. In our discussion, we will give a particular emphasis to the most signifi- 
cant differences between the CMSSM- and NUHM-seesaw scenarios. 

The numerical results presented in this section are mainly devoted to the particular case 
of Titanium nuclei, given that one expects a notable improvement of future experimental 
sensitivities in that case [13]. However, some additional predictions for other nuclei are also 
included here, for comparison. The case of Gold nuclei is of particular interest, since at 
present the most stringent bound is that of CR(// — e, Au) [10]. 

For the purpose of numerical evaluation, we begin by defining the input parameters at 
the gauge coupling unification scale, Mx- In the case of a CMSSM-seesaw scenario, and 
instead of scanning over the full (M1/2, Mq, Aq, tan/3, sign//) parameter space, we study 
specific points, each exhibiting distinct characteristics from the low-energy phenomenology 
point of view. We specify these parameters by means of the "Snowmass Points and Slopes" 
(SPS) [31] cases defined in Table 1. 



SPS 


Ml/2 (GeV) 


Mo (GeV) 


Ao (GeV) 


tan f3 


fi 


la 


250 


100 


-100 


10 


> 


lb 


400 


200 





30 


> 


2 


300 


1450 





10 


> 


3 


400 


90 





10 


> 


4 


300 


400 





50 


> 


5 


300 


150 


-1000 


5 


> 



Table 1: Values of M1/2, Mq, Aq, tan/3, and sign /x for the SPS points considered in the analysis. 



In the case of the NUHM-seesaw scenario, and in order to reduce the number of input 
parameters, we set Mq = M1/2 = Msusy, and explore the (Msusy,^o, tan/3, sign//, 5i, ^2) 
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parameter space considering the following intervals: 

250 GcV <MsusY < 1000 GeV, 
-500 GeV <Ao < 500 GcV, 
5 < tan 13 < 50, 
-2 < (5i , (52 < 2 . (31) 

In addition, we also consider the two possibilities, sign(/x) = ±1. 

To obtain the low-energy parameters of the model (and thus compute the relevant phys- 
ical masses and couplings), the full one-loop RGEs (including the neutrino and sneutrino 
sectors) are firstly run down from Mx to the right handed neutrino scale tjim- At this scale 
we impose the boundary condition of Eq. (9). After the decoupling of the heavy neutrinos 
and sneutrinos, the new RGEs are then run down from uim to the EW scale, at which the 
conversion rates are computed. Notice that, in the case of hierarchical heavy neutrinos, the 
sequential running is done from Mx down to mjvg and from m^Vi down to the EW scale. 
This implies that we do not take into account the running effects from the intermediate right 
handed neutrino scales, i.e. from m^Vs to m^Va and from mjvj to niNi- We have estimated 
these threshold effects by means of the leading logarithmic (LLog) approximation, verifying 
that they are indeed negligible for the numerical values chosen in the present work. 

The numerical implementation of the above procedure is achieved by means of the public 
Fortran code SPheno2 .2.2 [32]. The value of Mx is derived from the unification condition of 
the SU{2) and U{1) gauge couplings (systematically leading to a value of Mx very close to 
2 X 10^^ GeV throughout the numerical analysis), while l/x] is derived from the requirement 
of obtaining the correct radiative EW symmetry breaking. The code SPheno2 .2.2 has been 
adapted [33] in order to fully incorporate the right-handed neutrino (and sneutrino) sectors, 
as well as the full lepton fiavour structure. The computation of the /i — e conversion rates in 
nuclei, as well as of other LEV observables, has been implemented into the code by means 
of additional subroutines. 

Regarding the light neutrino masses and the C/mns matrix elements, we take the following 
input values: 

Am^„i = 8 X 10"^ eV% Am^^„^ = 2.5 x 10"^ eV\ m^^ = 10"^ eV, 

^12 = 30°, ^23 = 45°, ^i3<10°, (5 = = 02 = 0, (32) 

compatible with present experimental data (see, for instance, the analyses of [18]). We do 
not address the impact of non-vanishing C/mns phases (Dirac or Majorana) in the /i — e 
conversion rates. 
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Finally, although not used in this work, it is clarifying to recall that a simplified estimation 
of the generated fiavour mixing in the slepton sector can be obtained by means of the LLog 
approximation. Using the latter, the relevant off-diagonal slepton mass matrix element for 
the processes involving lepton flavour violation in the /i — e sector (as is the case oi /i — e 
conversion in nuclei) can be given as 

(Am|)2i = --^ (3 + Al) (yj L Y,),, ; L« = log ( ^) 5^ . (33) 

Writing {Y^LYy)2i using the parameterisation of Eqs. (5-10), and considering the limit of 
rriy^ = 0, 01 2 = 5 = (which is appropriate for the subsequent discussion), one obtains the 
following expression: 

vUy^ly^U = 

L33 rriNs [ci3 ^ym^ cl si S12 + ^/m^^ cl c\ S13] 

[y^Ci C2 Ci3 S23+ V^C2 Si (Ci2 C23 - S12 S13 S23)] 
+ L22 VTLj^.^ [^m^2 Cl3 (ct C3 - Si S2 S3) -512 - (C3 S\ + Cl S2 S3) S13] 

\-^/rn^^cvj, (ca si + ci S2 sa) S23 + ^fm^^ (ci C3 - si S2 S3) (ci2 C23 - S12 Si3 S23)] 

+ Lii mivi [- v^^cia {c\ s\ s\ + c\ S3) S12 + (-c^ c\ s\ + s* S3) S13] 

{^Jvrtjj^ Ci3 (-Cl C3 S2 + si S3) S23 - ^/rn^^ (cs si S2 + ci S3) (ci2 C23 - S12 S13 S23)] ■ 

(34) 

In what follows, we begin by investigating the theoretical predictions for the ji — e con- 
version rates in Titanium nuclei within the CMSSM-seesaw. 



4.1 Universality: CMSSM-seesaw 

The numerical results of the CR(/i — e, Ti) within the CMSSM-seesaw scenario are displayed 
in figures 2 through 5. The following discussion is focused on the most relevant parameters, 
namely mjVj, ^1,2,3, ^13, tan/?, Mq and Mi/2- 

In Fig. 2, we display the prediction of CR(// — e, Ti) as a function of the heavy neu- 
trino masses for the various SPS points, and for the particular choice di — {i — 1,2,3) 
and ^13 = 5°. We also consider the case of degenerate and hierarchical heavy neutrino 
spectra (respectively left and right panels). In both scenarios for degenerate and hierar- 
chical heavy neutrinos, we find a strong dependence on the the heavy neutrino masses. 
We also see that the rates for the various SPS points exhibit the following hierarchy. 
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BR4 > BRib > BRia > BR3 > BR2 > BR5. This behaviour can be understood in terms of 
the growth of the CRs with tan /3, and from the different mass spectra associated with each 
point. 

In the case of degenerate heavy neutrinos, we find the expected fast growing behaviour of 
CR(/x — e, Ti) as a function of the common neutrino mass uin- For the values oi uin within 
the studied interval [10^ GeV, 10^^ GeV] , the predictions for the CR(// — e, Ti) range over ten 
orders of magnitude. We also see that, for the chosen input parameter values, the predicted 
rates cross the experimental bound for the large uin region. In the latter, the Yukawa 
couplings can be large (for instance, Ygg and 1^32 can be while 1^22 ^21 can be of 

O(10~^)), leading to excessively large rates, so that these large niN values are disfavoured 
by data. The experimental bound is saturated for ttin values ranging from 2 x 10^^ GeV for 
SPS 4 up to about 10^^ GeV for SPS 5. In the case of hierarchical heavy neutrinos a similar 
behaviour of the predicted rates is found with respect to the heaviest neutrino mass, m^vg ■ 
We have also checked that the conversion rates do not significantly depend on m^Vi and m^Va, 
provided that their values are kept well below niNa- With the planned future sensitivity of 
10~^^ it will be possible to reach into wider regions of the heavy neutrino spectrum. Heavy 
neutrino masses above 10^^ GeV can be probed for the several considered scenarios. 

For most of the studied points, the previously illustrated dependence of the rates on the 
heavy neutrino masses is in agreement with the expected behaviour jm^v log mjvp obtained 
in the LLog approximation (as derived from Eq. (34)). However, a clear departure from this 
approximation is found for some points, the most remarkable being the case of SPS 5. This 
failure of the LLog approximation has been known to happen in some scenarios, for instance 
those with either large Aq, or low Mq and large M1/2 [34]. 

The predictions for CR(/i — e, Ti) as a function of the i?-matrix angles, ^1,2,3, are displayed 
in Fig. 3. In this case we have fixed the other relevant parameters as ^13 = 5°, ttin — 10^^ GeV 
and mjVj = (10^°, 10^^, 10^^) GeV (degenerate and hierarchical heavy neutrinos, respectively) 
and chosen SPS la. To fully explore the variation of the rates with the complex angles^ 9i, 
we have scanned the intervals < \9i\ < n rad and < arg^j < | rad. From this figure we 
see that the dependence on the three 9i is very similar in the degenerate case, whereas the 
same does not occur for hierarchical heavy neutrinos. In the former, the rates smoothly grow 

^Complex 6i may imply the presence of CP violation in the neutrino Yukawa couplings. In addition to 
affecting the LFV rates, these phases will induce contributions to flavour-conserving CP violating observ- 
ables, as is the case of charged lepton electric dipolc moments (EDMs). Throughout the present study we 
have verified that the associated predictions for the charged lepton EDMs are in agreement with current 
experimental bounds [2]. 
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Figure 2: CR(/i — e, Ti) as a function of the relevant heavy neutrino mass: m^r (on the left) and rriN^ (on the 
right), respectively associated with the degenerate and hierarchical cases. The predictions for SPS la (dots), 
lb (crosses), 2 (asterisks), 3 (triangles), 4 (circles) and 5 (times) are included. On the upper horizontal axis 
we display the associated value of {Y^)^^. In each case, we set 6*13 — 5°, and consider the limit where R ~ 1 
{9i = 0). A dashed (dotted) horizontal line denotes the present experimental bound (future sensitivity). 

with both modulus and argument, and are independent of 9i in the real case. In the latter, 
the rates are almost independent of 9^, and present a different minima pattern regarding 9i 
and ^2- The deep minima occuring in the real case are a consequence of the corresponding 
minima appearing in the relevant elements of the Yukawa couplings (as given by Eq. (9)). 
Notice that the observed behaviour of CR(/i — e, Ti) as a function of 9i can be indeed easily 
understood from the simple analytical expression obtained in the LLog approximation (cf. 
Eq. (34)). 

The most important outcome from Fig. 3 is that for both cases of degenerate and hierar- 
chical heavy neutrinos, complex values of 9.i can increase the /i — e conversion rates by almost 
five orders of magnitude with respect to the 6'j = case. Only for a few specific choices of 
9i (for instance real 9i or 6'2, in the hierarchical case) can we observe a strong decrease with 
respect to the 6'j = case, but clearly this is not a generic situation. 

In the following, and in order to simplify the analysis with respect to the other parameters, 
we will set 6^, = 0, and assume that the corresponding predictions for the CR(/i — e, Ti) will 
constitute a representative case for the lowest conversion rates. 

In Fig. 4 we show the dependence of the fi — e conversion rates on the light neutrino 
mixing angle 9is. The other parameters are set to ttln = 10^^ GeV, itln^ = (10^°, 10^^, 10^^) 
GeV, ( respectively for degenerate and hierarchical heavy neutrinos) and 9i = 0. All the SPS 
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Figure 3: From top to bottom, CR(/i — e, Ti) as a function of \9i\ [i = 1,2,3), for arg^i = 
{0, 7r/8, 7r/4, 37r/8, it/2] (dots, crosses, asterisks, triangles and circles, respectively). Both \6i\ and arg6'i 
are given in radians. On the left we consider degenerate heavy neutrinos (with rajq = 10^^ GeV), while on 
the right the hierarchical case is displayed (with ruNi = (10"'^°, 10"'^^, 10^'^) GeV). In all cases we take 6*13 = 5°, 
and set the CMSSM parameters to the SPS la case. A dashed (dotted) horizontal line denotes the present 
experimental bound (future sensitivity). 



points in Table 1 have been considered. For degenerate heavy neutrinos, the dependence 
on is softer than what is observed for the hierarchical case, leading to a variation in the 
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Figure 4: CR(^ — e, Ti) as a function of ^13 (in degrees), for SPS la (dots), lb (crosses), 2 (asterisks), 3 
(triangles), 4 (circles) and 5 (times). On the left we consider degenerate heavy neutrinos (with tun = 10^^ 
GeV), while on the right the hierarchical case is displayed (with niN. = (10^^,10^^,10^'*) GeV). In both 
cases we choose R = 1 {9i — 0). A dashed (dotted) horizontal line denotes the present experimental bound 
(future sensitivity). 



rates of at most one order of magnitude in the studied range of 0° < ^^13 < 10° (the only 
exception being SPS 5, where the variation can reach up to two orders of magnitude). In 
contrast, this figure clearly manifests the very strong sensitivity of the CR(/i — e, Ti) to the 
^13 mixing angle for hierarchical heavy neutrinos. In the hierarchical case, a variation of 613 
in the studied interval leads to an increase in the conversion rates by as much as five orders of 
magnitude. This huge variation is due to the strong decrease of this observable for very small 
^13 angles, as can be easily understood from the dependence on this angle of the dominant 
{LssiriNs s/mu^clclsis) term in Eq. (34). Furthermore, the minimum of CR(/x — e,Ti) is 
expected to occur at a vanishing mixing angle, but this being the value at the seesaw scale, 
i.e., 6'i3(mjv/) = 0. The deep minima in Fig. 4 are at Oi^{mz) — 0.2°, which is precisely the 
RGE shifted value at mz from 6i2,{mM) = 0. As 613 grows, the predictions for SPS 4, SPS 
la and SPS lb cross the present experimental bound. In particular, notice that for SPS 4, 
and for the present choice of input parameters, 6*13 values larger than 2° would be excluded 
by present data. 

An equally remarkable sensitivity to ^13 has been found in other fi — e violating processes, 
like /X — 67 and /i 3e, and also in tau decays as is the case of r ^ ej and r ^ 3e [34]. 
This interesting behaviour with ^13 has been proposed in [34] as a powerful tool to test the 
seesaw-I hypothesis for neutrino mass generation and, in case of a measurement of these 
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branching ratios, as a unique way to derive some hints on the seesaw parameters, especially 
on the value of m^vg. The fi — e conversion rates here presented will certainly add new 
interesting information on this type of analysis. Fig. 4 also shows that with the expected 
future sensitivity of 10~^^, the full 0° < ^^13 < 10° interval can be thoroughly covered. 

In the following study we will restrict ourselves to the hierarchical case where we have 
found this strong sensitivity to ^13. For definiteness, we will also fix the heavy neutrino 
masses and 6*13 to "reference" values of rriNi = (10^°) 10^^, 10^^) GeV and 6*13 = 5°. 
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Figure 5: Contributions to CR(/x — e, Ti): total (dots), 7-penguins (diamonds), Z-penguins (asterisks), 
-ff-penguins (crosses) and box diagrams (times). On the left we present the dependence on tan/3, for 
Mq — Ml/2 ~ 250 GeV and Aq = 0. On the right, we exhibit the evolution as a function of Mo{= ^^1/2)1 for 
tan /3 = 30 and ^0 = 0. In either case, we consider hierarchical heavy neutrinos with 111^. = (10^", 10^^, lO^"*) 
GeV, and set ^13 5°, and i? = 1 (6*^ = 0). A dashed (dotted) horizontal line denotes the present 
experimental bound (future sensitivity). 



Fig. 5 illustrates the predictions for the CR(/i — e, Ti) as a function of tan P, Mq and M1/2. 
Here we have separately displayed the various contributions to the fi — e conversion rates in 
order to conclude about their relative importance in this CMSSM-seesaw scenario. We set 
the values of the remaining CMSSM parameters to Mq = M1/2 = 250 GeV in the study with 
tan/3 (left panel) and to tan/? = 30 in the study with MgusY = Mq = M1/2 (right panel), 
taking = in both cases. We choose our "reference" values of ttin^^^ .^ = (10^°, 10^^, 10^"^) 
GeV, ^13 = 5°, and 9i = 0. 

In both panels of Fig. 5 we clearly observe the dominance of the photon-mediated contri- 
butions, which are in fact indistinguishable from the total CR, for all the explored parameter 
ranges. The dependence of the various contributions on tan/5 illustrates the expected fast 
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growing behaviour with tan^ (5 of the Higgs-mediated contributions, and the milder tan^ (5 
dependence of the photon-mediated ones. In addition, we see that the Z boson-mediated 
and the box diagram contributions are almost independent of tan/3. Although not displayed 
in this plot, we have also verified that the Higgs-mediated contribution is largely dominated 
by the exchange of H^, which is indeed the Higgs boson with enhanced couphngs to charged 
leptons in the large tan (5 regime. 

The decoupling behaviour for large Msusy of each of these contributions (CR-y, CRz, 
CKh and CRbox) is clearly manifested in the right panel of Fig. 5. The most important 
conclusion from this figure is that, within a CMSSM-seesaw scenario, the 7-penguin diagrams 
completely dominate the conversion rates, even for the largest tan/3 considered (tan/3 =50). 
Therefore, the total CR(// — e, Ti) does not manifest the Higgs contributions, so that in 
this universal scenario there is no chance for the /i — e conversion process to provide any 
information on the Higgs sector. We will see next that the situation is remarkably different 
in the non-universal case, where the Higgs contributions turn out to be much larger than in 
the universal case. 

4.2 Non- universality: NUHM- seesaw 

The numerical results for the NUHM-seesaw scenario are collected in figures 6 through 11. 

In order to study the influence of the hypothesis of non-universal Higgs soft SUSY break- 
ing masses, Mh^^, on the — e conversion rates, we have first explored the impact of the 
non-universality parameters 5i and 82 on the predicted Higgs boson masses. The values for 
these parameters have been taken to lie within the interval —2 < 81^2 < 2. 

The predictions for the relevant Higgs boson mass, m^/o, as a function of 5i and 82 are 
summarised in Fig. 6. We have chosen here the largest value of tan/3 = 50 and three 
representative values of Msusy = 250, 500 and 850 GeV for moderate, heavy and very 
heavy SUSY spectra, respectively. The other parameters are set to our "reference" values of 
mjv, = (10^°, 10l^ 10^^) GeV, Oi = 0, ^ = 0, ^13 = 5° and sign(/i) = +1. 

First, it is important to mention that not all the considered values of the 5i 2 parameters 
and Msusy allow for a correct SU{2) x U{1) breaking. In fact some particular choices for 
5i, (^2, and Msusy lead to unacceptable negative values of (and hence, negative m\o)- 
For instance, this is the case when 5i 2 are simultaneously positive or negative. Some other 
points, despite leading to a proper SU{2) x U{1) breaking, are nevertheless not acceptable, 
since they lead to a Higgs boson sector which is too hght, with masses below the present 
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Figure 6: Mass of the heaviest Higgs scalar (m^o) as a function of the non-universahty parameter ^i, for 
fixed values of 62 = {0, 0.5, 1, 1.5} (respectively dots, crosses, asterisks, triangles). The universahty case 
Si,2 = is represented by a large circle. We also take ruNi = (10^", 10^^, 10^^) GeV, set 9i = 0, Aq = 0, 
tan/3 = 50 and impose the relation Mq — Mi/2- The first three plots correspond to Mq = 250, 500 and 
850 GeV, respectively. On the fourth plot, we display the fi parameter as a function of the non- universality 
parameter 5i, for fixed values of 62, and for Mq — M1/2 ~ 500 GeV. 



experimental lower limits. To ensure that our results are indeed experimentally viable, we 
have included in this, and in the following figures, only the solutions where the three neutral 
Higgs boson masses are above the experimental bound for the lightest MSSM Higgs boson, 
which at present is 110 GeV for tan/5 > 5 (99.7% C.L.) [2]. The most interesting solutions 
with important phenomeno logical implications are found for negative 6i and positive 62, the 
choice selected for Fig. 6. In this figure, for all the explored values of 61 and we find a 
value of mffo that is significantly smaller than what one would encounter in the universal 
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Figure 7: On the left, mass of the heaviest Higgs scalar {171^0) as a function of the SUSY scale (Msusy = 
Mo = M1/2), for fixed values of Si = {—1.8, —1.6, —1, 0} (respectively crosses, asterisks, triangles and 
circles). We take toat, = (10^°, 10", 10^'^) GeV, and set e^ = 0, Aq = 0, tan/3 = 50 with 6*13 5°. On the 
right, the /x-parameter is displayed as a function of AfsusY, for the same choices of SUSY-seesaw parameters. 



case (here represented by the choice 5i = 82 = 0). This is truly remarkable in the case of 
large soft breaking masses, as can be seen, for instance, in the panel with Msusy = 850 
GeV, where low values of mfjo are still found, even close to the experimental limit. For 
completeness we have also shown in Fig. 6 the predictions for the fi parameter as a function 
of Si and 62- This parameter turns out to be nearly independent of 61, and its largest values 
are obtained for 62 = 0. 

The behaviour of the predicted niHO and /i parameter as a function of Msusy = Mq = 
Ml/2 is shown in Fig. 7. Here the specific values of 61 = {—1.8, —1.6, —1, 0} and ^2 = 
have been considered. This figure again illustrates the interesting departure from the linear 
behaviour of mno with Msusy, which is generic in the universal case {61,2 = 0). In contrast, 
the /i parameter conserves a similar linear behaviour with Msusy in all the studied scenarios 
(universal and non-universal). 

As a representative example of these interesting non-universal points, we explicitly refer 
to the case with 61 = —1.8 and 62 = 0, where the predicted masses are mno = 113, 174 and 
127 GeV for Msusy = 250, 500 and 850 GeV, respectively. For completeness, we have also 
collected the corresponding masses of the other relevant SUSY particles in Table 2. Notice 
that, in the case of Msusy = 850 GeV, this table illustrates a very heavy SUSY spectrum, 
even with a considerably heavy lightest SUSY particle, m^o = 362 GeV, but where the 
relevant Higgs boson is still light, 171^0 = 127 GeV. 
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MSSM masses 


MsusY (GeV) 


(GeV) 


250 


500 


850 


mr 


175 


415 


734 




258 


511 


867 


mr 


258 


511 


867 




307 


594 


985 




309 


607 


1025 


^l6 


323 


609 


1031 




281 


571 


971 


17102 


297 


601 


1022 




299 


605 


1028 




185 


395 


687 


^Xo 
Xl 


379 


679 


1075 


Xi 


99 


207 


362 


X2 


185 


394 


687 


X3 


363 


668 


1067 


Xi 


377 


678 


1074 


rriho 


110 


119 


123 




113 


174 


127 



Table 2: Relevant MSSM spectra for Mq = M1/2 = Msusy, tan/3 = 50, Aq = 0, Oi = 0, $13 = 5°, 
MjVi = (10^°, lO^S 10") GeV, 61 = -1.8 and 62 = 0. 



In the following we present the predictions of the /i — e conversion rates in Titanium 
nuclei within the NUHM-seeesaw scenario. First we display in Fig. 8 the CIi{fi — e, Ti) as a 
function of Mq = M1/2 = Msusy and of Aq for the particular choice 61 = —1.8 and 62 = 0. 
In order to illustrate the impact of the non-universality hypothesis on the conversion rates, 
we have separately displayed in this plot the various contributions from the 7-, Z-, Higgs 
mediated penguins and box diagrams. We observe a very distinct behaviour with Msusy 
of the Higgs-mediated contributions when compared to what was found for the CMSSM 
(universal) case, shown in Fig. 5. In fact, for the choice of input parameters in this plot, 
the Higgs-mediated contribution can equal, or even exceed that of the photon, dominating 
the total conversion rate in the large Msusy region. Both photon- and Higgs-mediated 
contributions are similar around Msusy = 700 GeV. These larger Higgs contributions are 
the obvious consequence of the lighter Higgs boson mass values encountered in this region, as 
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previously illustrated in Figs. 6 and 7. The non-decoupling behaviour of the SUSY particles 
for the large MgusY regime can be seen in the Higgs contribution, and thus in the total rates 
for the Higgs-dominated case. 

For completeness, we have also explored other choices of Aq and sign(/i). The case 
of sign(yu) = —1, whose numerical results are not presented here, does not evidence any 
interesting new feature. In fact, there is a much more reduced 6i, 62 parameter space allowing 
for the correct SU{2) x U{1) breaking. In addition, for sign(/i) = —1 we have not found 
solutions displaying as small values of m/^o as in the case of sign(/i) = +1. The predicted 
Higgs contributions to the conversion rates are correspondingly smaller, and therefore less 
interesting. Regarding Aq, the right panel in Fig. 8 shows that all the contributions are 
essentially independent of the value of the universal trilinear coupling, so that our selected 
value, Aq = 0, is in fact a good representative point. 
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Figure 8: Contributions to CR(/i — e, Ti) as a function of Mo{— M1/2) (left) and Aq (right): total (dots), 7- 
penguins ( diamonds), Z-penguins (asterisks), iJ-penguins (crosses) and box diagrams (times), for 5i — —1.8 
and 82 = 0. We set tan/3 = 50 and take 6*13 = 5°, i? = 1 [9, = 0) and rriN, = , 10", lO^'') GeV. On the 
left = 0, while on the right we choose Mq{= M1/2) = 700 GeV. In each case, a dashed (dotted) horizontal 
line denotes the present experimental bound (future sensitivity). 



Within the NUHM-seesaw scenario, we have also studied the ji — e conversion rates for 
other nuclei. The case of Gold is particularly interesting since its present experimental bound 
of 7 X 10^^^ [10] is more stringent than the present bound for Titanium (4.3 x 10^^^ [9]). In 
Fig. 9 we display the predicted /i — e conversion rates for Al, Ti, Sr, Sb, Au and Pb, as a 
function of MgusY- We have chosen two light, two moderate and two heavy nuclei and we 
have fixed the other parameters to those of the previously elected non-universality reference 
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point (with 6i = —1.8 and ^2 = 0). For completeness, the values of the relevant parameters 
for these nuclei, Z^s, Fp and Fcapture, have been collected in Table 3 and follow [35]. In 
this figure we clearly see that throughout most of the explored Msusy interval, the relative 
conversion rates obey the hierarchy CR(/i — e, Sb) > CR(/i — e, Sr) > CR(yU — e, Ti) > 
CR(/i — e, Au) > CR(/i — e, Pb) > CR(/i — e, Al), in agreement with the generic results 
in [35]. We do not find a significant difference in the large Msusy region, where the Higgs- 
contribution dominates the ratios. The predicted rates for Ti, Au and Pb tend to converge 
whereas the corresponding curve for Al nuclei deviates slightly from the others at large 
-^susY, but we do not consider these differences among the predictions for the various nuclei 
to be relevant. The most important conclusion from Fig. 9 concerns the fact that we have 
found predictions for Gold nuclei which, for the input parameters in this plot, are clearly 
above its present experimental bound throughout the explored Msusy interval. However, it 
should be recalled that the formulae here used for these estimates come from approximations 
that may not properly work for the case of very heavy nuclei. These heavy nuclei deserve a 
more dedicated and refined study. 
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Figure 9: ^ — e conversion rates for various nuclei as a function of Mq = M1/2 in the NUHM-seesaw. We 
display the theoretical predictions for Sb, Sr, Ti, Au, Pb and Al nuclei (diamonds, triangles, dots, asterisks, 
times and crosses, respectively). We have taken m^Vi — (10^", 10^^, 10-^'*) GeV, Aq = 0, tan/3 = 50, ^13 = 5° 
and R—\ {di ~ 0). The non-universality parameters are set to Si ~ —1.8 and 62 ~ 0. From top to bottom, 
the horizontal dashed lines denote the present experimental bounds for CR(/x — e, Ti) and CR(/i — e, Au). 
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^Nucleus 






rcapt(GeV) 


?^A1 


11.5 


0.64 


4.64079 X 10" 


19 


iTi 


17.6 


0.54 


1.70422 X 10- 


18 


ISr 


25.0 


0.39 


4.61842 X 10- 


18 




29.0 


0.32 


6.71711 X 10- 


18 


]fAu 


33.5 


0.16 


8.59868 X 10- 


18 




34.0 


0.15 


8.84868 X 10- 


18 



Table 3: Values of Z^ff, Fp and Fcapt for different nuclei, as taken from [35]. 

Before proceeding with our analysis, let us briefly mention that for the region investigated 
in Fig. 9, the SUSY contributions to the anomalous magnetic moment of the muon, = 
range from a^^^^ = 10~^ for Msusy = 250 GeV to a^^^^ = 10"^, in association with 
-^susY = 850 GeV. The latter values are in fair agreement with the observed excess in 
when compared to the SM prediction, which, at the 3.8 a is given by a^^^^ = a^^^ — — 
3.32 X 10~^ at 3.8(7 (for a review, see for instance [36] and references therein). 



To complete our study of the ji — e conversion rates in the NUHM-seesaw scenario wc have 
compared the theoretical predictions for the CR(/i — e, Ti) with those for the BR(/i 67). 
We recall that both observables are sensitive to the same leptonic mixing given by the slepton 
mass matrix entries connecting the first and the second generation. In the usual photon- 
penguin dominated case, the latter two quantities are known to be highly correlated, and 
this is indeed what occurred for the CMSSM-seesaw discussed in Sec. 4.1. In other seesaw 
scenarios, as for instance, SUSY-GUT seesaw [37] or the inverse seesaw [38], this strong 
correlation still persists. However, for some scenarios where the photon-mediated diagrams 
are no longer the dominant contributions to the conversion rates, the strong correlation 
between CR(// — e, Ti) and BR(// — > 67) can be lost. For instance, this loss of correlation 
has been found in the case of Littlest Higgs Models, as recently pointed out in [39] . 

We have also found an interesting loss of correlation in the present case of the NUHM- 
scenario, where, as previously discussed, the Higgs-contributions can be the dominant ones. 
The departure from the strongly correlated regime for (CR(/i — e, Ti), BR(/i — > 67)) 
is illustrated in Fig. 10, considering several choices of the neutrino mixing angle ^13 = 
10°, 5°, 1°,0.2°. For aU plots the predictions for (CR(// - e, Ti), BR(// 67)) have been 
derived for several choices of the non- universality parameter 5i, scanning over the following 
interval 250 GeV < Msusy < 1000 GeV. In each of the panels, the predictions for (CR(// — e. 
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Ti), BR(/i — > 67)) that correspond to Si — S2 — fall upon a straight hne, which strongly 
supports the correlated behaviour of the two observables in this case. As Msusy increases 
within the considered interval, (CR(/x — e, Ti), BR(// 67)) moves left and downwards 
along the straight hne due to the obvious decrease of the rate with Msusy- 

However a clear departure from the previous strongly correlated predictions is found for 
other values of 5i, 52- In particular, for the specific 5i and 82 values where, as previously 
shown, the Higgs contributions dominate the /i — e conversion rates, the predicted (CR(/x — e, 
Ti), BR(/x — > 67)) points exhibit a different behaviour, deviating from the straight line 
associated with the universal case. The separation between the correlated and uncorrelated 
regimes is maximal for the Si — —1.8, ^2 = non-universal case, as can be clearly understood 
from our previous results. We find this loss of correlation a very promising phenomenon that 
could be fully explored if future sensitivities of 10~^^ are reached. 

Secondly, it is clear from Fig. 10 that even in the most pessimistic situation of very small 
^13, the theoretical predictions for CR(/i — e, Ti), and in particular the corresponding curved 
line, are well above the horizontal hne at 10~^^. This is quite a challenging possibihty, since 
for those high values of Msusy ~ 850 GeV, whose predictions he at the left end of the curved 
and straight lines, the predicted BR(// — > 67) is far below the planned 10~^^ sensitivity. This 
clearly reflects that /i — e in nuclei can be a very competitive process to study LFV within 
the SUSY seesaw. 

Finally, and to summarise the most striking results for these NUHM-seesaw scenarios, 
we plot in Fig. 11 the ratio of the two predicted rates, CR(/x — e, Ti)/BR(/x 67) as a 
function of m^jo. Since both observables exhibit the same dependence on mjVj, di and ^13, 
the consideration of this ratio of rates allows to reduce the number of relevant parameters 
to tan/9, Msusy and 61^2- These last two are clearly the leading ones given that they drive 
the solutions to the interesting low ttiho values. In this figure, and in order to maximise the 
Higgs-contribution to the total /i — e conversion rates we have again considered the extreme 
tan P — 50 value. For consistency, we have set the remaining parameters to their reference 
values, but as we have said, they will not play a relevant role in this study. 

Leading to this scatter plot, we have scanned in the intervals — 2 < (5i < 0, < ^2 < 2 
and 250 GeV < Msusy < 1000 GeV. The most important conclusion from this plot is that 
the ratio CR(// — e, Ti)/BR(// 67) can deviate from the constant prediction of 5 x 10~^ 
of the universahty case by as much as a factor of almost 10. For the scan here conducted, 
the maximum value of this ratio of rates is found for 5i — —1.7, 82 — 0.1 and Msusy = 876 
GeV, and its size is 0.04. 
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Figure 10: CR{fi - e, Ti) versus BR(^ ^ ej) for 250 GeV < AfguSY < 1000 GeV, and Si = 
— 1.8, —1.7, —1.6, (crosses, triangles, asterisks, dots, respectively). We set S2 = 0, and take rriNi — 
(10^°, 10", lO^") GeV, Ao = 0, tan/3 = 50 and i? = 1 {e^ = 0). From left to right and top to bottom, the 
panels are associated with 6*13 — 10°, 5°, 1° and 0.2°. In each case, the horizontal and vertical dashed (dot- 
ted) lines denote the present experimental bounds (future sensitivities) for CR(/i — e, Ti) and BR(// — > 67), 
respectively. 

Considering larger values of MgusY and identical intervals for 5i 2 leads to somewhat 
similar results: one finds the same pattern of clusters departing from the constant value of 
the universal case, but the maximum value of CR(/i— e, Ti) /BR(/i —>■ 67) is in general smaller 
than the 0.04 obtained in the scan of Fig. 11. The reason why this ratio is not improved at 
larger values of Msusy than 1 TeV is because the acceptable solutions producing the proper 
SU{2) X U{1) breaking do not lead to sufficiently light Higgs bosons. 

Even without the knowledge of the seesaw parameters, a measurement of CR(/i — e, Ti) 
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Figure 11: Ratio CR(^ — e, Ti)/BR(/i 67) as a function of the Higgs mass, niffo. We take niNi = 
(10^°, 10", 1014) GeV, Ao = 0, tan/3 = 50, 6*13 = 5° and i? = 1 {9^ = 0), and scan over 250 GeV < 
-^SUSY < 1000 GeV, —2<Si< 0, and < 62 < 2 (grey dots). We have highlighted specific choices of 
82 — 0, 0.1, 0.2, 1.1 (crosses, triangles, diamonds, dots, respectively). In each tilted cluster, we have also 
indicated the values of 61 associated with the 62 coloured points. The universality limit {5i = 82 = 0) is 
denoted by a circle. Asterisks denote points with Afgugy — 876 GeV and 62 — O-l. 



and BR(/i —>■ 67), together with information on tan/5 and the SUSY scale, may allow to 
shed some light into the Higgs sector. 



5 Conclusions 

In this work we have extensively studied the lepton flavour violating process of yU — e con- 
version in nuclei, within the context of the SUSY-seesaw. In particular, we considered two 
distinct scenarios, the CMSSM-seesaw, and the NUHM-seesaw, obtained by partially re- 
laxing the universality conditions of the Higgs boson masses. Throughout our analysis, we 
compared our theoretical predictions with the present experimental bounds, and with the 
challenging future sensitivities. In fact, the latter may convert processes like CR(/i — e, Ti) 
into one of the most sensitive probes to new physics. 

We have presented here the first full one-loop computation of the fi — e conversion in 
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nuclei, including the complete set of SUSY-loop diagrams: 7-mediated, Z- and Higgs-boson 
mediated penguins and box diagrams. We have also provided the full analytical results 
working in terms of physical eigenstates (for all intervening SUSY and Higgs particles) . 

For the CMSSM-seesaw, we have considered the dependence of the conversion rates on 
the several parameters defining the scenario. Choosing the well known SPS benchmark 
points to specify the CMSSM parameters, we focused on the most relevant parameters in 
the neutrino sector, namely on the heavy neutrino masses (m^vj, the complex Qi mixing 
angles and the still undetermined angle of the C/mns matrix, ^13. As discussed here, the CRs 
exhibit a very pronounced dependence on the previous parameters, with variations that can 
reach up to ten orders of magnitude in the case of mjv^ and up to five orders of magnitude in 
the cases of Qi and ^13, for the investigated ranges. In turn, this strong dependence implies 
that a comparison of the theoretical predictions with the present experimental bound allows 
to derive indirect upper bounds for the unknown seesaw parameters. 

We have pointed out that the highest sensitivity is found for the case of hierarchical 
heavy neutrinos. In this case, the conversion rates are essentially dependent on mjvg and 
^1^2, manifesting an extreme sensitivity to ^13 (for the case of vanishing Q^. In fact, the values 
of these parameters in the upper part of their studied intervals, 10^^ GeV < m^^ < 10^^ GeV, 
< 1^1,2! < TT, < arg(^i_2) < 7r/2 and 0° < 613 < 10° are already in conflict with the present 
upper bounds on CR(// — e, Ti) and CR(// — e, Au). 

We have put special emphasis on the sensitivity of the CR(/i — e,Ti) to ^13, given that 
either a measurement, or a more stringent bound on this parameter is expected in the near 
future [40]. Therefore, and once ^13 is measured, a dedicated study of the — e conversion 
rates could provide some insight into the potentially unreachable heavy neutrino parameters. 

In all the studied examples of the CMSSM-seesaw, the dominant contribution to the 
II — e conversion rates clearly arises from the photon-penguins. Even though we have verified 
that the Higgs contributions do indeed grow with tan^ P, they induce contributions which 
are several orders of magnitude below those of the photon (which grow as tan^ /3) for all the 
studied interval 5 < tan/3 < 50. A very interesting departure from this situation occurs 
when one relaxes the universality condition for the Higgs soft breaking masses, and this 
fuelled our interest to consider the NUHM-seesaw. 

In the case of the NUHM-seesaw, we explored the influence of the non-universality hy- 
pothesis of the soft SUSY breaking masses 2 on the /i — e conversion rates. The Si and 
S2 parameters that describe the departure from universality in the Higgs sector have an im- 
portant impact on the predicted Higgs boson masses. In particular, we have found regimes 
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for 5i^2 with very interesting phenomenological implications, namely the possibility of a light 
Higgs spectrum, even in the limit of large soft SUSY masses. As a concrete example, we 
recall that for the reference choice of 5i — —1.8, S2 — 0, we find uiho — 113, 174 and 127 
GeV for Msusy = 250, 500, 850 GeV respectively (in turn associated with moderate, heavy 
and very heavy sparticle spectra). 

The distinctive NUHM-seesaw scenarios associated with light bosons and a relatively 
heavy SUSY spectra induce very interesting and unique predictions for the /i — e conversion 
rates. Specifically, we have shown that in the large Msusy = Mq — M1/2 region (e.g. 
above 700 GeV), there is a strong enhancement in the Higgs-dominated rates, leading to a 
remarkable loss of correlation between the CRs in nuclei and the BRs of // — > 67 decays. 
As we aimed at illustrating in Fig. 10, the departure from the linear correlation of these 
two observables can be sizable. It is worth stressing that if both these rates and 613 are 
measured, values of BR(// — > 67) and CR(// — e, Ti) that clearly deviate from the expected 
SUSY-seesaw ratio in the photon-dominated case, can provide indirect information into the 
structure of the Higgs sector. 

It is also important to remark that with the expected future sensitivities, /i — e conversion 
in nuclei maybe sensitive to LFV signals that lie beyond the reach of the future sensitivities 
to /X — > 67 decays. For example, this can occur for a heavy SUSY spectrum, and very small 
values of ^13. 

Finally, we considered the predictions for the ratio CR(/x — e, Ti)/BR(/i — 67) as a 
a function of ttiho in NUHM-seesaw scenarios, comparing the results with those obtained 
for the CMSSM-seesaw case. The most important conclusion to be drawn from this study 
(which is presented in Fig. 11) is that in the NUHM-seesaw one can observe a clear deviation 
from the constant prediction of the CMSSM-seesaw by as much as a factor close to 10. If 
such deviation is indeed observed, we can obtain some indirect hints regarding the SUSY 
Higgs sector. 

In summary, with the expected future sensitivities, /i — e conversion in nuclei can clearly 
be more competitive for the study of LFV in SUSY-seesaw than // — > 67, and certainly 
provide an important tool for the study of the Higgs sector. 
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A One-loop formulae for /i — e conversion in nuclei 

In this appendix wc collect all the analytical results of the SUSY one-loop diagrams that 
contribute to the — e conversion rates in nuclei. These are summarised by the photon-, Z- 
boson- and Higgs-boson- penguins and box diagrams in Figs. 12, 13, 14 and 15, respectively. 
In the following subsections we present the relevant formulae for each separate contribution. 
All the loop functions in the formulae are taken from [33,41]. 

A.l Form factors for the 7//e vertex 




7 7 7 

Figure 12: Relevant SUSY one- loop diagrams for the photon- mediated contributions to /i — e conversion in 
nuclei. 

Our convention for the form factors defining the 7//e vertex is as follows: 

ie [g^a(Af Pi + A^Pr) + im^a^0^{A^^PL + A^Pr)] , (35) 
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where q is the off-shell photon momentum, P^^r — (1 ^75)72, e is the electromagnetic charge 



and is the muon mass. 



In the SUSY-seesaw context there are one-loop contributions to these form factors that 
come from the chargino and neutralino sectors respectively, 



j^L,R ^ ^{n)L.R ^ a = 1, 2 . 

The neutralino contributions are given by, 
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where xax = 'mio /mf and the indices are ^4 = 1, .., 4, X = 1, .., 6. 
The chargino contributions are given by 
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where in this case xax — m^_/m|^ and the indices are A — 1,2, X — 1,2,3. Notice that 
in both neutralino and chargino contributions a summation over the indices A and X is 
understood. 
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A. 2 Form factors for the Zjie vertex 

Our convention for the form factors F^ r defining the Z^e vertex is as follows: 

-^7M [FlPl + FrPr] . 



(43) 




Figure 13: Relevant SUSY one- loop diagrams for the Z-mediated contributions to /U — e conversion in nuclei 



The Z-boson form factors have also the two kinds of contributions, from neutralinos (n) 
and charginos (c), 

FLiR) = Fi^) + F^l^^ . (44) 
The results for the corresponding form factors are the following: 
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where again the indices are A,B — 1,..,4, X,Y — 1,..,6 in the contributions from the 
neutralino sector and A,B — 1,2, X,Y — 1,2,3 in the contributions from the chargino 
sector, and a summation over the various indices is understood. 
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A. 3 Form factors for the Hjie vertex 

Our convention for the form factors r defining the H^e vertex is as follows: 

[HlPl + HrPr] . 



(49) 




Figure 14: Relevant SUSY one-loop diagrams for the Higgs-mediated contributions to /x — e conversion in 
nuclei. 



As in the previous cases, we separate the contributions from the neutralino and chargino 
sectors, 
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The results for the form factors are the following, 
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= H 



-mlBi{mlo^,ml)N^AxN;^Xx + ^/.^x^5o(m|o , mf^)iV2i;,7V^^;^ 

mem^5i(m|o^, my iVf^^iV^^;^ + mem^o^5o(m|o^, m?^)iV,ViV5:^] } , (51) 

(52) 



ip) 

L,n 



p=l,2,3. 



Correspondingly, the result for the chargino contribution H^^j^^ ^ can be obtained from the 



previous H^j^^ ^ by replacing everywhere. 



jyL{R) 

Dl{r) 



u 

X~ 

qL{R) 

Wl{r) 



In the previous formulae, the index p refers to the each of the Higgs bosons. Concretely, 
Hp — hP,H^,A^ for p — 1,2,3, respectively. The other indices are again A,B — 1,..,4, 
X,Y — 1, .., 6 in the contributions from the neutralino sector and A,B — 1,2 and X,Y — 
1, 2, 3 in the contributions from the chargino sector. A summation over all the indices is also 
understood. 
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A. 4 Contributions from box diagrams 



We follow here a simmilar notation as in the previous formulae for the separate contributions 
from the neutralino and the chargino sectors. Our convention for the box diagrams at the 
quark level is iBg. 




Figure 15: Box diagrams contributing to /x — e conversion in nuclei. 



The results for the vector contributions are the following: 
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(53) 
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The results for the scalar contributions are given by: 
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(58) 
(59) 



(60) 
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The indices in the previous formulae are again, A,B = 1, ..,4, X,Y = 1, ..,6 in the contri- 
butions from the neutralino sector and A,B = 1,2, X,Y = 1, 2, 3 in the contributions from 
the chargino sector. A summation over all the indices is also understood. 



B Relevant couplings for fi — e conversion in nuclei 

In this appendix we collect the formulae for the couplings that are relevant in this work. 
We follow the same notation for the couphngs as in [33] and include here some of the 
formulae presented there, for completeness. The couplings are expressed in the physical 
eigenstate basis, for all the MSSM sectors involved: sleptons Ix {X — 1, ..,6), sneutrinos ux 
{X — 1,2,3), neutralinos Xa = charginos Xa — 1>2) and the neutral Higgs 

bosons Hp {p — 1,2,3) — hP,H^,A^. Notice that the case Hp — is given for completeness 
but, as explained in the text, it does not contribute to the /i — e conversion in nuclei in the 
coherent approximation assumed in the present work. 

The notation for the SM parameters that appear in the following couplings is as follows: 
g is the SU{2) gauge coupling, m/ is the fermion mass, mwi t^z are the VF-boson and 
Z-boson masses, respectively, and 6w is the weak angle. 
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B.l Neutralino couplings 



The couplings for neutralinos that enter in the one-loop diagrams computed here are the 
following: 

= -9V2 1 2^^;:^3 ^ ^l3^[?,3,5)x + tan 9wN*^,Rfl,,^x } . (63) 

= { 2mrcos /? ^^3^('-3,5)x + I tan gv^7Vl,<?,,e)x } . (65) 

TV^W = -^^^{-^ (-^tang^TV^i + TV^,^ i^;?;3,)^ + ^^^cos/j ^-^^^S'^^)^} '^^^^ 
TVfH ^ -^v^l^^^^^AT^^i?!^^^^^^^^ (67) 

Here, arc the 6x6 rotation matrices for the charged slcpton, down squark and 

up squark sectors, respectively, and is the 4x4 rotation matrix for the neutralino sector. 
For completeness, we have written the full set of couplings, including the three fermion 
generations. The displayed notation for the sfermion rotation matrices with three entries 
R{ , ,) correspond with the three generic possibilities to fermion index i. The fermion masses 
are correspondingly, m^. = rUe, rrin, mr] nidi — ^d, nis, mi, and m„. = m„, mc,mt. Notice also 
that, although we use the same notation for the squark and slepton sectors, and since we 
have not included mixing in the quark sector, the 6x6 rotation matrices R^'^^ and i?(")are 
block diagonal in flavour space and only L — R mixing occurs in that case. 
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B.2 Chargino couplings 



The couplings for charginos that are present in the one-loop diagrams computed here are 
the following: 
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(70) 
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(74) 



Here i?*^^^ is the 3x3 rotation matrix for the sneutrino sector, and U and V arc the 2x2 
rotation matrices in the chargino sector. The displayed notation for the three entries in the 
sfermion rotation matrices is as in the previous neutralino couplings. The rotation matrices 
for neutralinos and charginos can be found in [42] and [43]. 



B.3 Z boson couplings 
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Zvx^ coupling: 



Zll coupling: 



Zqq coupling: 
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We have used here the short notation sw — sin^^^ and cw — cos 9w- 



B.4 Higgs boson couplings 

HpX\x% couphng: 
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with 
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with A; = {AiY"- / {YiY'' (at the EW scale), i = 1, 2, 3 for Z = e, /x, r, respectively. 
Hpi'x^Y coupling: 



^XY — 9 9lL,v^\X -"■ir "I" 9lL,v^2X ^2Y "I" 9lL,u^3X ^3Y 



with 



coupling: 



(p) 



*M pM I ^(p) p*('^) 0(1^) 



(p) 
9ll.v 



St] = 'J 



2 cos 



c(p) 



2mvK COS (3 



(p)* 
1 ) 



Hpdd coupling: 



?(p) 



'^''^ ^2mM/cos/?'^^ 



(p)* 



c(p) _ c 
'^R,d — '-' 



(p)* 



ifpMM coupling: 



S 



(p) 



c(p) 



-9 



2mvK sin /3 



:(J' 



(p)* 
2 ) 



r,(p)* 



43 



In all the above equations, 
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We have also used here the standard notation for the low-energy MSSM soft-gaugino-mass 
parameters Mi 2 and the /i parameter. 
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